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We study the rotational predissociation of atom - molecule complexes with very small 
binding energy. Such complexes can be produced by Feshbach resonance association 
of ultracold molecules with ultracold atoms. Numerical calculations of the predis- 
sociation lifetimes based on the computation of the energy dependence of the scat- 
tering matrix elements become inaccurate when the binding energy is smaller than 
the energy width of the predissociating state. We derive expressions that represent 
accurately the predissociation lifetimes in terms of the real and imaginary parts of 
the scattering length and effective range for molecules in an excited rotational state. 
Our results show that the predissociation lifetimes are the longest when the binding 
energy is positive, i.e. when the predissociating state is just above the excited state 
threshold. 



I. INTRODUCTION 



The energy levels of ultracold atoms can be shifted by an external magnetic field. This 
gives rise to shifts of bound energy levels of atom - atom complexes. Under certain conditions, 
an adiabatic variation of the external field can be used to move a scattering state of ultracold 
atoms to below threshold. This results in the formation of atom - atom complexes with 
extremely small binding energy (< 10" 8 eV). This technique - known as Feshbach-resonance 
association - is widely used for the production of ultracold molecules from ultracold atomsi~-. 
Feshbach resonance association can also be used to link ultracold atoms with ultracold 
molecules or ultracold molecules with ultracold molecules^ to produce weakly bound atom 

- molecule or molecule - molecule complexes. The binding energy of such complexes can 
be tuned from zero to about 0.5 cm" 1 by varying an external field. If the molecules are 
initially in a ro-vibrationally excited state, these weakly bound complexes may undergo ro- 
vibrational predissociation. While the ro-vibrational predissociation of van der Waals atom 

- molecule complexes has been the subject of numerous studies^—, little is known about 
the predissociation of atom - molecule complexes with binding energies < 0.5 cm -1 . In 
the present paper, we explore the dependence of the rotational predissociation lifetimes of 
weakly bound atom - molecule complexes on the binding energy in the limit of zero binding 
energy. 

The present work is motivated by two recent research developments. First, the cre- 
ation of ultracold molecules in the ro-vibrational ground state, achieved by several research 
groups^i^, makes the measurements of near-threshold rotational predissociation possible. 
The predictions of the present work can therefore be readily tested in experiments with 
ultracold molecules. Second, several recent studies showed that an ensemble of ultracold 
molecules trapped on an optical lattice may undergo collective rotational excitations anal- 
ogous to Frenkel excitons in solid-state molecular crystals^"—. Ultracold molecules on an 
optical lattice can therefore be used as a model for quantum simulation of energy transfer 
in solid-state crystals and molecular aggregates^. Of particular interest would be a study 
of energy transfer in the presence of a disorder potential. A disorder potential can be 
introduced to an ensemble of molecules on an optical lattice by adding atoms to certain 
lattice sites. If the atoms are coupled to molecules by Feshbach resonance association, the 
atom - molecule interactions may introduce controllable perturbations affecting the Frenkel 
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excitons. The rotational predissociation of atom - molecule complexes would, however, lead 
to loss of molecules from the optical lattice, thus limiting the experiments. It is therefore im- 
portant to explore if the lifetimes for the rotational predissociation can exceed the timescales 
for the rotational energy transfer between molecules on an optical lattice^. 

The dependence of the predissociation lifetimes on the binding energy can be qualita- 
tively understood using perturbation theory^— . Using the unperturbed single channel 
radial functions for the bound and continuum states, it can be shown that the predissoci- 
ation lifetime is inversely proportional to the square root of the binding energy. However, 
perturbative calculations fail to give quantitative results due to the difficulty of representing 
the near-threshold bound state wavefunctions and the neglect of couplings between contin- 
uum states. An alternative approach is based on multichannel effective range theory and 
the calculation of elastic and inelastic scattering matrix elements for the excited state in 
the limit of zero collision energy^^. This method can be used to estimate the location 
and lifetime of the most weakly bound predissociating state. However, all previous cal- 
culations were based on the assumption that the couplings between molecular states are 
exceedingly small, resulting in exceedingly weak inelastic scattering. In the vicinity of a 
Feshbach resonance, giving rise to a predissociating state near threshold, this assumption is 
generally not valid. In this paper, we extend the work of Balakrishnan and coworkers^ and 
Forrey and coworkers^ to derive the relations between the predissociation lifetime and the 
scattering parameters of the atom - molecule complex without any assumptions about the 
strength of inelastic scattering. We calculate the predissociation lifetimes for the prototyp- 
ical Mg( 1 S , )+NH( 3 E _ , n = 1) ultracold system, where n is the molecular rotational angular 
momentum. We find that near threshold the predissociation lifetime can be as large as 100 
fis. 



II. THEORY 

At low collision energies, the scattering cross section of an atom and a molecule is entirely 
determined by the s-wave contribution to the partial wave expansion and the scattering 
observables can be expressed in terms of a few parameters. For simplicity of discussion, we 
assume that the atom is structureless. Using effective-range theory^, the first few terms in 
a power series expansion of the cotangent of the scattering phase shift Si, or equivalently the 
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diagonal S- matrix element Su = exp(2iSi), for a molecular state i, can be written as 



where fc, is the scattering channel wave number defined by h 2 kf = 2fi(E — Ei), [i is the 
reduced mass, E is the total energy and Ei is the energy of the molecular state %. Each 
molecular state gives rise to a scattering channel. The constants Oj and T{ are the channel 
scattering length and effective range, respectively. The atom-molecule interaction potential 
depends asymptotically on the atom - molecule separation R as R~ n . The scattering length 
exists if n > 3 and the effective range exists if n > 5. 

If state i is the absolute ground state of the molecule, only elastic scattering is possible, 
\Sa\ = 1 and both the scattering length and the effective range are real. If state i is an 
excited state, inelastic scattering becomes possible, \Su\ < 1 and both the scattering length 
and the effective range become complex—^—, a% = oti — ifii and = tr^ + ir\^. Indeed for 
two s-wave states, a ground and an excited state, the energy dependence of the scattering 
observables near the excited state threshold can be characterized by just six real parameters: 
the real and imaginary parts of the scattering length and effective range of the excited state, 
the tangent of the ground state scattering phase shift at threshold and an additional effective 
range parameter that describes the first order energy dependence of the ground state phase 
shift at threshold. 

A Feshbach resonance occurs when a bound state corresponding to an asymptotically 
closed channel > E) is the same as the collision energy of lower-energy molecular states. 
In the single channel the system is tuned across a resonance, the S"-matrix element 

describes a unit circle in the complex plane and the scattering length exhibits a pole when 
Su passes through —1. In the multichannel case, Su no longer describes a unit circle in the 
complex plane and thus the scattering length no longer exhibits a pole but an oscillation^. 
The location of the Feshbach resonance corresponds to a pole in Su. Regarding Su as a 
function of the complex energy the pole occurs at 



where E& is the location of the bound state/resonance and T is the resonance width, related 
to the bound state lifetime by r = h/T. 

To allow analysis below threshold (at negative channel energies), we define the complex 
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8 = E B - i- 
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channel wave number Ki = —iki such that 

Balakrishnan et al.— related the location and lifetime of the bound state/resonance below 
threshold to the complex scattering length calculated above threshold. In subsequent work^, 
the treatment was expanded to include the complex effective range. However in this latter 
work only the limiting case of small inelasticity was considered, in which q>^ and r^j ^> 
rij. In the vicinity of a Feshbach resonance, the pole in the real part of the scattering length 
is suppressed^ 1 ^ and thus the magnitude of /3 may be comparable with or greater than that 
of a (see Figure 1, for example). 

The location of a Feshbach resonance also corresponds to a pole in the scattering ampli- 
tude f{k). Following Balakrishnan et al.— , the scattering amplitude can be written as 



/(*) 



1 ~i -i 

1 T i .2 , 



(4) 



di 2 

for a potential that asymptotically decays faster than R~ 4 . In the limit of Ki — > 0, the pole 
occurs when Ki = 1/dj and thus when 

2fia 

and the bound state energy and width are explicitly 



r /i|a| 4 ' 1 ' 

In the single-channel case ((3 = 0), as the bound state goes from just below to just above 
threshold, at which point the bound state becomes virtual, the scattering length goes through 
a pole changing sign from +oo to — oo^. In the presence of inelastic channels ((3 ^ 0) the 
bound state is coupled to the other channels which leads to a shift in the bound state energy. 
As the bound state is tuned from just below threshold to just above, the bound state crosses 
threshold (E-q = 0) when a 2 = (3 2 and becomes a virtual state (r < 0) when a < (or 
< 0). In the presence of inelastic channels the true bound state can thus exist at energies 
greater than the asymptotic threshold energy. 

Incorporating the effective range as a correction to the pole in f(k) occurs when 



Making no assumptions about the nature of the scattering, the complex energy at which the 
pole occurs is 

Introducing the notation 

p + ia = (r R + in) (a + i/3) = ar R - 0n + i(an + /3r R ) (10) 

and 
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the complex energy at which the pole occurs can be explicitly written as 

e = _l_ K - r!) - 2ir R n ( 1 _P+^_ SR + iS \ 
H |rj| 4 V \ a \ J 

where the location and width of the bound state are 
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III. NUMERICAL CALCULATIONS 

To test the validity of equations d5J) and (jUj), we numerically calculate the location and 
width of the highest energy s-wave bound states of a Mg( 1 S') atom bound to a rotation- 
ally excited NH( 3 S~) molecule. The potential energy surface of the Mg+NH complex was 
calculated by Soldan, Zuchowski and Hutson^. This potential energy surface was used in 
quantum mechanical scattering calculations^ demonstrating that the sympathetic cooling 
of NH by Mg has a good prospect of success. The Mg-NH energy surface is anisotropic, 
providing significant couplings between the predissociating and continuum states. 

The rotational structure of NH exhibits fine structure due to the coupling of the molecular 
spin s to the rotational angular momentum n. The molecular energy levels are characterized 
by the quantum number of the total angular momentum j = n + s of the molecule. There is 
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thus a single ground rotational state \n = 0, j = 1) and three states \n = 1, j = 0, 1, 2) cor- 
responding to n — 1. The end-over-end rotation of the complex is described by the angular 
momentum L which couples to j to form the total angular momentum J = j + L. The ma- 
trix of the total Hamiltonian describing the Mg - NH complex can be represented in the fully 
coupled angular momentum basis set \(ns)jLJM), as described in^i. We perform the scat- 
tering calculations using the MOLSCAT package^, in which sets of coupled equations are 
constructed for each J and parity (— l) n+L+1 in the fully coupled basis set \(ns)jLJM) with 
n < 6 and L < 12. The coupled differential equations are integrated using the hybrid log- 
derivative method of Alexander and Manolopoulos^ based on a fixed-step-size log-derivative 
propagator in the short-range region (2.5 to 50 A) and a variable-step-size Airy propagator 
in the long-range region (50 to 250 A). The log-derivative solution is matched to the asymp- 
totic boundary conditions^ to obtain the S-matrix. To calculate the scattering length and 
effective range for use in equations © and <Q, the S'-matrix is calculated above threshold 
over the energy range 10 pK to 100 nK and then fitted to equation (CQ). 

To vary the location of bound states with respect to threshold, we introduce a scaling 
factor to the potential 



In the following, we consider the first (most weakly bound) predissociating state below the 
Mg+NH(ri = l, j — 1),L — threshold. Figure [TJ shows the complex scattering length 
and effective range as a function of A sca i e . As the bound state energy crosses threshold 
(at A sca i c ~ 1.03), the real and imaginary parts of the scattering length show the expected 
oscillation and peak respectively. The real part of the effective range also shows a suppressed 
pole as a passes non-resonantly through zero. An additional resonance is also present at 
Ascaie ~ 1-065, due to an n = 2 state crossing threshold. 

The location and width of each predissociating bound state are determined by calculating 
the S'-matrix as a function of energy. For a single open channel, the scattering phase shift 
Si goes through n at the bound state energy following the Breit-Wigner energy dependence 



where Sbg(E) is a slowly varying background term. In the multichannel case, the single 
channel phase shift no longer follows to Breit-Wigner form and it is the S-matrix eigenphase 
sum a(E)&, the sum of the phases of the eigenvalues of the S-matrix, that follows the Breit- 
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scaled 



(R) = A scak y(R). 
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FIG. 1. The real and imaginary parts of the Mg+NH(ra = l,j = 1) s-wave scattering length 
a = a — if} and effective range r = tr + in as functions of the scaled interaction potential, are 
shown as a (red, solid line), f} (blue, dashed line), (green, dashed-dotted line), and r\ (purple, 
dotted line). 



Wigner form. When the bound state exists well below threshold, V < E t h IC s — E-q, the width 
and location are found by fitting the eigenphase sum across resonance to the Breit- Wigner 
equation using the RESFIT computer program^, as shown for A sca i c = 1.1 in figure |2j 
When the bound state exists close to threshold T > |i?thres — Eb\ or just above threshold 
Eb > -Ethrcs, the threshold cusp prevents the eigenphase sum from going through the full 
Breit- Wigner form. An example is shown in figure [2] for A sca i = 1.035 at which point the 
predissociating state is just below threshold. The location and width of the resonance must 
thus be extrapolated from the eigenphase sum below threshold. Although the calculated 
eigenphase sum values are essentially exact, this limits the accuracy of the locations and 
widths obtained from the close coupling calculations for bound states very near (below or 
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E(cm" 1 ) E(cm" 1 ) 

FIG. 2. Left panel: The location and width of a predissociating bound state at energy E-q well below 
the threshold energy E^\ iVCS is calculated by fitting the Breit-Wigner form (blue, solid line) to the S- 
matrix eigenphase sum (EPS), obtained from solving a set of coupled-channel (CC) equations (red, 
crosses). Right panel: For a predissociating state close to (r > Eb — £thres) or above threshold, as 
shown for A sca i c = 1.035, the location and width of the state are extrapolated from the eigenphase 
sum below threshold. 

above) threshold. 

The binding energy, the width and the lifetime of the most weakly bound state are shown 
in Figure |3] as functions of the real part of the scattering length. It can be seen from the 
lower panel of Figure [3] that the calculation based on equation ([9]) agrees with the close 
coupling numerical calculations far better than the calculation based on equation (jSJ) for 
a > 50 A when the binding energy is negative. As the bound state moves above threshold 
the extrapolation of and T from the eigenphase sum below threshold breaks down and 
the calculated results become inaccurate. The gap in the data is due to the presence of the 
additional n = 2 Feshbach resonance at A sca i e ~ 1.065. 

Starting from the bound state existing below threshold, as we decrease A sca i e , increasing 
a, the binding energy decreases and the predissociation lifetime increases. When a ~ (3 « 
160 A, the bound state crosses threshold, which is represented by the sharp increase of the 
curves in the upper panel of Figure El Decreasing A sca i e further, decreases a and moves 
the bound state above threshold. At low scattering lengths, the width of the state shows 
a dramatic decrease before becoming negative, at which point the bound state becomes 
virtual. At these small positive scattering lengths, it can be seen that the Mg-NH bound 
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state lifetime can be on the order of tens to hundreds of microseconds. From Eq (7) T — > 
as a — > + , however with the effective range correction Eq (14), r — > when a ~ 18 A. 
Thus in principle the width of the state can be arbitrarily small. 

The results presented in Figure 2 are specific for the Mg - NH complex. However, Eqs. 
(13) and (14) are universal. To explore the universality of the predissoication lifetimes 
presented in Figure 2, we repeated the calculations for different values of the interaction 
potential scaling parameter A sca i e . As A sca i e increases, the strength of the interaction potential 
increases. This increases the number of bound energy levels supported by the potential. 
The scattering length as a function of A sca i c exhibits a series of peaks corresponding to zero- 
energy bound states (see the middle panel of Figure 3). The upper panel of Figure 3 displays 
the predissociation lifetimes as functions of the scattering lengths calculated using different 
ranges of the scaling parameter corresponding to different peaks in the middle panel of Figure 
3. The different curves in the upper panel of Figure 3 thus correspond to different interaction 
strengths. It can be seen that that the predissociation lifetimes of the bound states very 
near threshold are sensitive to the interaction strength. This is because the magnitude of 
the effective range is different for different interaction strengths, as demonstrated in the 
lower panel of Figure 3. However, the predissociation lifetimes remain similar at small 
values of the scattering length. In particular, the large values of the predissociation lifetimes 
corresponding to the bound state above threshold appear to be insensitive to the interaction 
strength. This indicates that the large magnitude of the predissociation lifetimes just above 
the excited state threshold is a general phenomenon. 

The presence of an atom weakly bound to a molecule perturbs the rotational energy 
structure of the molecule. Atom - molecule Feshbach resonances can be used to introduce 
controllable perturbations to the rotational structure of ultracold molecules. This could be 
used to create an ordered array of ultracold molecules on an optical lattice in the presence 
of tunable disorder potential. It is therefore important to quantify the effect of the atom 
- molecule potential on the rotational structure of the molecules in Feshbach-resonance 
associated atom - molecule complexes. 

The energy separation of the NH(n = 0, j = 1) rotational ground state and NH(n = 
l,j = 1) can be approximated as 

2A 

(E nlil - Koji) « 2B mt n(n + 1) - 7sn + (17) 
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where in this approximation the rotational constant is B rot = 16.34685 cm" 1 rather than the 
spectroscopic value of 16.343 cm -1 and the spin-rotation 7 sn and spin-spin A ss constants are 
-0.055 and 0.920 cm -1 respectively^. When the Mg atom is present the energy separation 
is modified, resulting in an effective rotational constant B* ot . Figure shows the energy 
difference between B* ot and B TOt calculated using Eq. OH]), as a function of a. It can be seen 
that the rotational constant of the molecule can be effectively modified by linking the atom 
and the molecule with a Feshbach resonance. 

IV. SUMMARY 

We have shown that multichannel effective range theory can be used to calculate the 
location and lifetime of the most weakly bound predissociating state near threshold, even 
in the presence of strong couplings between molecular states (strong inelasticity). As the 
predissociating state is tuned closer to and moves across threshold, the lifetime of the state 
increases. We have found that the predissociation lifetime of the states just above the 
excited state threshold is the largest and can be tens to hundreds of microseconds long. 
For illustrative purposes, the location of the bound state was tuned in the present work 
by varying the depth of the potential. In experiments, bound states can be tuned across 
threshold using externally applied electric and magnetic fields^. 

An ensemble of molecules on an optical lattice can be used to study collective excita- 
tions characteristic of solid-state crystals^"—. The presence of an atom weakly bound to 
a molecule creates a perturbation in the rotational energy structure of the molecule. Our 
calculation shows that this perturbation can be tuned from zero to ~ 100 MHz. Adding ul- 
tracold atoms to certain lattice sites and using Feshbach resonances could thus be exploited 
to modify the rotational structure of selected molecules trapped on an optical lattice. Ro- 
tational excitation of molecules in such a system would resemble collective excitations in a 
crystal with impurities^. The characteristic time for transfer of rotational excitation energy 
between molecules in different sites of an optical lattice is ~ 10 /is. The calculations of the 
present work show that the lifetimes for the rotational predissociation can be as large as ~ 1 
ms. 

Mixtures of ultracold atoms and molecules have recently been created in several experi- 
mental laboratories^ 1 ^. The ultracold atom - molecule complexes undergoing the rotational 
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predissociation studied in this work can be created by first linking ultracold molecules in 
the ground rotational state with ultracold atoms using a Feshbach resonance and then ex- 
citing the molecules to a rotationally excited state. The external field used for the Feshbach 
resonance association can be varied to adjust the binding energy of the atom - molecule 
complex. We hope that the results of the present work will stimulate the measurements of 
near-threshold rotational predissociation. 
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FIG. 3. The binding energy (upper panel) and the lifetime (lower panel) of the most weakly bound 
predissociating state as functions of the real part of the scattering length. The close coupling 
calculations are represented by diamonds, the results of equation ([5]) - by the blue, dashed line, 
and the the results of equation © - by the rec-^solid line. 




FIG. 4. Upper panel: The lifetime of the near-threshold predissociating state as a function of 
the scattering length modified by varying the interactional potential using different ranges of the 
scaling parameters corresponding to different peaks in the middle panel. Middle panel: The real 
part of the Mg+NH(n = 1, j = 1) s-wave scattering length as a function of the interaction potential 
scaling parameter. Lower panel: The real part of the effective range for the Mg+NH(n = 1, j = 1) 
collision complex as a function of the scattering length corresponding to different peaks in the 
middle panel. 
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